Vortex-Peierls States in Optical Lattices 
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We show that vortices, induced in cold atom superfluids in optical lattices, may order in a novel 
vortex-Peierls ground state. In such a state vortices do not form a simple lattice but arrange them- 
selves in clusters, within which the vortices are partially delocalized, tunneling between classically 
degenerate configurations. We demonstrate that this exotic quantum many-body state is selected 
by an order-from-disorder mechanism for a special combination of the vortex filling and lattice 
geometry that has a macroscopic number of classically degenerate ground states. 



The existence of quantized vortices is one of the most 
dramatic manifestations of the macroscopic wavefunction 
("order parameter") in superfluid (SF) Boseeases. Con- 
siderable theoretical [ij and experimental effort has 
thus been applied to the study of SF vortices in cold 
atoms systems. Vortices are usually regarded as classical 
objects, that form regular lattices due to long-range re- 
pulsive interactions between them. Intuitively, this point 
of view seems to be in accord with the fact that vortices 
are topological objects, i.e. the existence of an isolated 
vortex in the system can be established by observing 
phase winding infinitely far away from the vortex core. 
However, analogies with two-dimensional (2D) electronic 
systems in high magnetic fields, exhibiting the fractional 
quantum Hall (FQH) effect, suggest that under certain 
conditions vortex lattices may be melted by quantum 
fluctuations, and various strongly correlated vortex liq- 
uid states may thus emerge 3]. 

Optical lattices offer additional opportunities to ex- 
plore the quantum mechanical behavior of vortices by 
allowing one to tunc the strength of quantum fluctua- 
tions. Several approaches to stabilizing FHQ states in 
cold atoms systems using optical lattices have already 
been proposed 0, 0, 01 ■ More generally, of interest are 
situations in which vortices may behave as strongly in- 
teracting quantum particles, moving in a periodic opti- 
cal lattice potential. Such a possibility appears to be 
rather counterintuitive, since naively we think of vortices 
as macroscopic objects. Nevertheless, the common view 
nowadays is that vortices in 2D systems can be consid- 
ered as quantum particles with a finite mass It is 
then interesting to find experimentally observable phe- 
nomena in cold atoms systems in which this quantum 
nature of vortices is manifested. In particular, manifesta- 
tions of the quantum mechanical behavior of vortices may 
be most dramatic in situations when long-range intervor- 
tex interactions are frustrated, which strongly enhances 
the effect of quantum fluctuations. Such a frustration in 
optical lattice systems may be engineered by choosing the 
appropriate combination of the optical lattice geometry 
and the vortex density. We should point out that some 
classical commensuration effects between vortex lattices 
and the underlying optical lattice pinning potential have 
already been studied, both theoretically \p\ and experi- 



mentally [9j, but the possible quantum effects have not 
been previously considered. 

In this Letter we discuss an example of such a quantum 
mechanical behavior of SF vortices in an optical lattice. 
We show that for a particular combination of the optical 
lattice geometry and the vortex flUing, for which classical 
vortex configurations are strongly frustrated, a vortex- 
Peierls (VP) state is realized. By VP state we mean a 
vortex lattice, in which vortices are not localized at the 
maxima of the optical lattice potential, but are instead 
partially delocalized, resonating quantum mechanically 
between degenerate pinned configurations. Peierls order- 
ing has been extensively studied, most recently in the 
context of quantum magnetism of localized spin systems 
(see Ref. jlflj for review) . Here we demonstrate that such 
ordering may under certain conditions occur for vortices 
in optical lattice superfluids. 

We consider a SF system of cold bosonic atoms, loaded 
in an optical periodic potential with the dice lattice ge- 
ometry, shown in Fig^ The fascinating features of the 
quantum mechanics of particles on the dice lattice in a 
perpendicular magnetic field were first pointed out by 
Vidal et al. and extensively studied in a number 

of subsequent works (see e.g.J2j and references therein). 
Atoms in an optical lattice can be described by the Bose- 
Hubbard model ll^ , in which bosons are assumed to tun- 
nel between nearest-neighbor sites of the lattice and in- 
teract when they are on the same site. We will assume 
that the average number of bosonic atoms per site of the 
dice lattice n is an integer. In this case increasing the 
on-site interaction energy U relative to the hopping am- 
plitude t until U jnt ~ 1 will induce a SF-Mott Insulator 
(MI) transition [l^ when n bosons will be localized on 
each site to minimize the interaction energy. 

To induce vortices in the SF we add an effective per- 
pendicular "magnetic field" . We will comment on par- 
ticular methods, that could be used to create such an 
effective field, at the end of the paper. We will focus on 
a specific value of the flux per plaquette of the dice lattice 
2nf = 27r/3, where / = 1/3 has the physical meaning of 
the vortex filling, i.e. the average number of vortices per 
plaquette. Centers of dice plaquettes can be associated 
with sites of the dual kagome lattice, which correspond 
to the maxima of the optical lattice potential. It is then 
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FIG. 1: (Color online). Optical dice lattice created by su- 
perimposing three laser field potentials /i(r), /2(r) and /3(r) 
(see text), ai = |a; — and a2 = |x + are the ba- 
sis directions. Lighter areas correspond to potential minima, 
biequivalent sites in the three-site unit cell of the dice lattice 
are labeled as 1,2 and 3. 

convenient to assume that the vortex cores are located 
on the kagome lattice sites. Since vortices interact via 
a long-range repulsive potential, they will try to arrange 
themselves in patterns on the kagome lattice, that maxi- 
mize the distance between each vortex and its neighbors. 
As shown by Korshunov at filling factor / = 1/3 
the set of vortex configurations that minimize the classi- 
cal interaction energy between the vortices consists of all 
states, where every triangular plaquette of the kagome 
lattice is occupied by exactly one vortex. The number of 
such configurations grows exponentially with the system 
size and the classical ground state of vortices has only al- 
gebraic order at zero temperature, but no true long-range 
order [l5| . 

Nevertheless, below we will demonstrate that strong 
quantum fluctuations near the SF-MI transition lift the 
classical degeneracy and select a vortex state with a true 
long-range order. This state has, however, a manifestly 
quantum mechanical nature, in that the vortices are not 
simply localized on sites of the kagome lattice, but are in- 
stead partially delocahzed over plaquettes, thus exhibit- 
ing what we call a VP ordering. 

To approach the problem analytically, we will take ad- 
vantage of the presumed proximity of our system to the 
SF-MI critical point. Since the bosons are at an integer 
filling, the standard Landau- Ginzburg (LG) theory can 
be applied not too far away from the critical point. The 
imaginary time LG functional of the SF order parameter 
fields is given by [lfij |: 



S ^ dr 
Jo 



+ Y,i\^r'^>^\' + rm' + uim^ + u,\<^>,f)]. (1) 

i 

Here is the local SF order parameter, /j, denote the 
nearest-neighbor vectors of the dice lattice, J ^ Uht, 
tt4,M6 > and r tunes the system across the SF- 
MI transition. To find the ordering patterns near the 



transition, we need to diagonalize the first term in 
Eq-O ,1Z]- Choosing Landau gauge for the vector po- 
tential Aifj, = iirfij-il — 5^x) one obtains the follow- 
ing dispersion for the lowest Hofstadter band: e(k) = 
— J [6 -I- cos(fci) -I- cos(fc2) — 2 cos(A:i — + •\/3sin(fci) 

1 /2 

+ -\/3sin(fc2)] , where k = fcibi + fc2b2 and bi_2 = 
ia; ± are the reciprocal lattice vectors of the dice 
lattice. The boson dispersion has two minima inside 
the first Brillouin zone of the dice lattice, at wavevec- 
tors ko = (0,27r/3) and ki = (27r/3,0). The corre- 
sponding eigenvectors are given by = (1,1,0) and 
— (e"^'^'/'^, 0, 1). We can then write the lattice order 
parameter fields as linear combinations of these two 
low energy modes $o-(ri) — X^^^o i '■Pf.'^G^^^'^'^' j where the 
index i labels the unit cells of the dice lattice, a = 1,2,3 
labels sites within each unit cell and ipi are the fields, 
corresponding to the low energy boson modes. To obtain 
the LG action in terms of the fields fg in its most general 
form we need to know how these fields transform under 
the symmetry operations of the dice lattice. The relevant 
operations are the elementary translations along the basis 
directions Ti and T2, rotations by 7r/3 around the six- fold 
coordinated sites i?7r/3 and reflections with respect to the 
X- and y-axes Ix,y These transformations are given by: 

T ■ ,0 ^ ,0 „-27rif/3 rp , ,n c27ri(^-l)/3 p 

li . (fi ^ (fiee ' , I2 ■ —>■ 'Pee ^ " , it^/s . 

,r, ,n „~2iri(e+l)/3 J ■ ,n . ,_*„-27ri€/3 r . , 

ft 'Pe+ie ^ " , ix ■ tpe ^ f^e , -ly ■ ^ 

V'l+ie^'^'^'^, where the subscripts of the fields are taken 
modulo 2. Using these transformations, the most general 
form of the imaginary time LG action is found to be ■ 
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In the mean-field approximation, which we expect to 
be accurate for the effective classical 2-|-l-dimensional 
system, described by Eq.©, three different SF phases are 
possible, depending on the signs of the v and w couplings. 



1. V 



> 0: either |i^o| 7^ or \ipi \ ^ 0. 



2. w < 0,w < 0: |(/5o| = 7^ 0- The relative phase 
B ~ arg[(pg(p]^] is determined by the last term in 
Eq.l^l) and is given by 9 = 27rn/3, n = 0, 1, 2. 

3. u < 0,w > 0: |<y5o| = Ifil 7^ 0. The relative phase 
is given by 6* (2n + l)7r/3, n = 0, 1, 2. 

To reinterpret the states we have found in the vortex 
language, it is convenient to calculate gauge-invariant 
supercurrents on each bond, which we define as J^^ ~ 
Im (<&,*<i>,;+^e~'^'''). In the v > state we find that su- 
percurrents vanish on every bond. This fact, combined 



FIG. 2: (Color online). (Left) Order parameter configuration 
corresponding to the v > state. Direction of each arrow 
represents the phase of the corresponding order parameter 
field, while length represents the magnitude. Expectation 
values of the order parameter vanish on all sites of type 3 
(shown by dots). (Right) Corresponding vortex configuration. 
Shaded triangles contain a vortex, uniformly delocalized over 
the three sites of each triangle. Note that the supercurrents 
vanish on every bond of the dice lattice. 




FIG. 3; (Color online). (Left) Order parameter configuration 
corresponding to the VjW < state. (Right) Corresponding 
supercurrent configuration. Plaquettes, containing vortices 
(shown by circles), are the ones that have supercurrents cir- 
culating around them in the counterclockwise direction. 

with the picture of this state in terms of the SF order 
parameter, leads to the vortex configuration shown in 
the right panel of Fig[21 Since this state contains con- 
figurations, in which two vortices are located on nearest- 
neighbor sites, it can not be the true ground state of the 
vortices. 

Calculating supercurrents in the second candidate 
ground state, realized when v,w < 0^ we obtain the con- 
figuration shown in the right panel of Fig|31 Vortices 
in this case are localized on the dice lattice plaquettes, 
which have supercurrents circulating around them in the 
counterclockwise direction. This configuration is a mem- 
ber of the classical ground state manifold, since none of 
the vortices have nearest neighbors. 

Finally, the supercurrent pattern in the state, realized 
when V < and w > 0, is shown in the right panel 
of FigEl One can see that none of the dice plaquettes 
in this case have a full vortex localized in it. Instead 
the vortices appear to be bound in partially delocalized 
triplets, populating the six dice lattice plaquettes, ad- 




FIG. 4: (Color online). (Left) Order parameter configuration 
corresponding to the n < 0, > state. Expectation value of 
the order parameter vanishes on some of the six-fold coordi- 
nated sites (shown by dots). (Right) Corresponding supercur- 
rent configuration. Note that none of the plaquettes contain 
a full vortex. Instead, vortices are bound into triplets pop- 
ulating the groups of six plaquettes, adjacent to six-fold co- 
ordinated sites with zero order parameter expectation value. 
Below, SF order parameter phase snapshots, corresponding 
to the two resonating vortex triplet configurations are shown. 

jacent to the six-fold coordinated sites with zero order 
parameter expectation values. The fact that the order 
parameter vanishes on these sites means that vortices are 
moving around such sites, locally destroying phase coher- 
ence. It also means that the boson number does not fluc- 
tuate on such sites, i.e. the bosons on these sites are in 
the Mott phase. The only vortex state that is consistent 
with this picture, and also does not violate the constraint 
of no-nearest-neighbor vortices is the one in which the 
vortex triplets resonate between two degenerate config- 
urations (corresponding instantaneous order parameter 
phase configurations are explicitly shown in Fig^J. This 
state is clearly more energetically favorable than the state 
in Figl^l since it does not violate the no-nearest-neighbor 
constraint but also allows the vortices to gain kinetic en- 
ergy by partially delocalizing over hexagonal plaquettes 
of the dual kagome lattice. We can estimate the energy 
gain in this state due to the vortex delocalization, com- 
pared to the w < state, as follows. Interaction energy 
of the vortices (energy of the supercurrents) is of order 
fit per lattice site. Vortex kinetic energy (energy gain 
from phase fluctuations) is of order U per site. There- 
fore, the energy gain in the w > state, compared to the 
w < state, is of order nt{U/nt)^, which is not small 
when U ~ fit. Thus, we find that the state in Fig0]is the 
ground state of the vortices in our problem. Since the 
vortices are not localized on sites of the kagome lattice, 
but are partially delocalized over plaquettes, we call this 
state a VP state. 

Let us now discuss how to observe the VP state ex- 
perimentally. Creating an optical lattice with a dice 
geometry experimentally is more difficult than most 
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other 2D lattices, but fortunately still possible with the 
current technology. We propose the following proce- 
dure (see Ref.J^i for an alternative proposal). One 
first creates a kagome lattice, which is done using 
the laser field potential proposed in Ref.^9|: /i(r) — 
Eti [cos (k, • r + ^) + 2 cos (ik, • r + ^)]\ where 
ki = (tt, \/37r), k2 = (tt, -VStt), kg = (-27r, 0), CTi = 
(73 = 1 and a2 = —1. Here the maxima of the potential 
correspond to the kagome lattice sites. To create a per- 
fect dice lattice we superimpose two additional laser po- 
tentials, that have a triangular lattice geometry: /2(r) — 
4ELi cos^ ■ and hir) = -4^^! cos^ (g, • r), 
where gi = (7r,7r/V3), g2 = (0, -27r/%/3), and gs = 
(— TT, tt/VS). Here /2(r) has maxima at the six-fold coor- 
dinated sites of the dice lattice, while I3 (r) has minima 
at both three-fold and six-fold coordinated sites. The su- 
perposition of Ii (r), /2(r) and /3(r) creates a perfect dice 
lattice, in which all potential wells have equal depth. Ef- 
fective magnetic flux in this setup can be created by the 
rotating mask method |^ , which generates a rotating op- 
tical lattice potential. 

Alternatively, the proposal of Ref. 4J| can be used to 
create the effective perpendicular field. One uses the 
combination of a time-dependent quadrupolar potential 
V{t) — Vqp 8iTi{ujt)xy and a temporal modulation of the 
tunneling amplitudes for different nearest-neighbor di- 
rections in the dice lattice. Modulation of the tunneling 
amplitudes can be achieved in our case by varying the 
strengths of the three components of the /3(r) poten- 
tial. To achieve the flux of 27r/3 per plaquette, the pa- 
rameters of the quadrupolar potential have to be chosen 
such that Vqp/hu) — In this setup, time-of-flight 

interference imaging can be used to study the periodic 
structure of the vortex lattice. Lattice periodicity, how- 
ever, is not enough to distinguish the VP state we found 
from the state in Fig|3| since they have identical recip- 
rocal lattice vectors. VP ordering can be detected by 
analyzing nontrivial noise correlations |20l |. that will be 
present in the time-of-flight image. Namely, in addition 
to the sharp peaks in the density distribution (i.e. flrst 
order correlation function), corresponding to the recipro- 
cal lattice vectors gi/3 of the votex lattice, one should ob- 
serve strong incoherent background. Such a background 
would be absent in a regular superfluid, in which all vor- 
tices are strictly localized. This is due to the fact that 
the vortex motion locally destroys phase coherence on 
some sites of the optical lattice. In the case of the VP 
state, described above, the incoherent background con- 
tains correlations at the same wavevectors gi/3, which 
may be found by measuring the density autocorrelation 
function. If the perpendicular field is created by the ro- 
tating mask technique, one can use interference between 
two identical co-rotating condensates jl^l to image the 
vortex lattice configurations. 



In conclusion, we have proposed that vortices in op- 
tical lattice SF may exist in VP ground states, which 
are direct analogs of valence-bond-solid states of inter- 
acting bosons. In particular, we have demonstrated that 
in the case of a dice optical lattice with vortex filling of 
1/3 per plaquette, the ground state of the vortices is a 
plaquette VP state, in which vortices bind into triplets, 
that resonate between two degenerate configurations on 
plaquettes of the dual kagome lattice. Such unconven- 
tional vortex ordering is a result of an order-by-disorder 
phenomenon, where extensive degeneracy of frustrated 
classical vortex configurations is lifted by quantum fiuc- 
tuations. 
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